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Abstract: Certain time-like singularities are shown to be resolved already in classical 
General Relativity once one passes from particle probes to scalar waves. The time evolution 
can be defined uniquely and some general conditions for that are formulated. The Reissner- 
Nordstrom singularity allows for communication through the singularity and can be termed 
"beam splitter" since the transmission probability of a suitably prepared high energy wave 
packet is 25%. The high frequency dependence of the cross section is However, 
smooth geometries arbitrarily close to the singular one require a finite amount of negative 
energy matter. The negative-mass Schwarzschild has a qualitatively different resolution 
interpreted to be fully reflecting. These 4d results are similar to the 2d black hole and are 
generalized to an arbitrary dimension d > 4. 
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1. Introduction 



The importance of singularities in General Relativity and their resolution is well appre- 
ciated. The most intriguing ones are considered to be those inside black holes, and the 
Big-Bang-like cosmological singularities. String theory offers some cases where a time-like 
singularity gets resolved: the orbifold, the flop, and the conifold. These are all singulari- 
ties in compact factors of spacetime, the last two occurring in Calabi-Yau manifolds, and 
their resolution often involves adding light matter ("twisted sectors," "wrapped D-branes") 
which lives on the singularity (for a review, see and references therein). Recently, there 
was a renewed effort to include time-dependence in string theory (for a review, see |2| and 
references therein). Ultimately, one would like to understand space-like singularities such 
as Schwarzschild and cosmological singularities, but it would be fair to say that there were 
no breakthroughs yet. 

The study described here was motivated by bold attempts to cross the singularity in a 
stringy black-hole model. The 2d black hole is a case where an algebraic coset construction 
requires one to add to spacetime additional manifolds which are glued together over the 
singularity (see figures ||,^) . For the uncharged 2d black hole it was found Q that waves sent 
from the "additional" infinity are fully reflected, while a charged 2d black hole does have a 
non-zero transmission amplitude across the singularity and, moreover, wave functions are 
smooth there [^. 

Here we show that there is nothing intrinsically stringy in those calculations ^. Ac- 
tually this interesting phenomenon happens once one considers the wave equation on this 
background, namely, by passing from point-particle probes ^ to waves. More precisely, even 
though the wave equation is singular, the ordinary differential equation (ODE) obtained 
after separating time, and whose solutions we denote by 4)^^, allows for a unique, smooth 
continuation through the singularity. Interestingly, this property is not special to the 2d 
charged black hole, but it generalizes to the 4d charged Reissner-Nordstrom (RN) black 
hole, as well as to d > 4 RN black-holes (moreover, for the 2d and 4d cases the relevant 
ODE is actually completely regular). This observation strongly suggests that the above 
classical RN spacetimes should be viewed as being made of two parts which are naturally 
connected across the singularity at r = 0. The two parts are (i) r > which is a usual 
(positive-mass) RN black hole, and (ii) r < 0, the region beyond the singularity, which is 
a negative-mass spacetime. 

Another motivation for considering the negative-r part of RN comes from the real- 
istic astrophysical black holes, which are known to be rotating in general [||, ^ ^. In 
the analytically-extended Kerr geometry (describing a stationary spinning black hole) , the 
Boyer-Lindquist coordinate r goes smoothly from positive to negative values. The region 
r < is an asymptotically-flat universe (distinct from the original positive-r external uni- 
verse). There is a curvature singularity at r = 0, but this singularity is a ring rather than 

^ There are a' corrections in string theory, which were computed in the 2d case. These give rise to extra 
phases in scattering amplitudes which vanish in the semi-classical limit. 

^The usual definition of a singular spacetime is "geodesic-incompleteness," namely, incomplete time- 
evolution for particle probes. 
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a hypersurface. A timelike geodesic heading towards r = will generically avoid the ring 
singularity, and smoothly pass "through the ring" into the negative-r asymptotic region. 
Obviously, wave packets will also pass through the ring, though with partial reflection. 
Spherically-symmetric charged black holes were often considered as useful toy models for 
the more realistic spinning black holes, due to the remarkable similarity in the inner struc- 
ture of the two black-hole types. From this perspective we may regard the RN solution 
analyzed here as such a simplified toy model for a spinning black hole. 

Initially one would expect the time evolution of a wave to be ill-defined in the presence 
of a singularity since some boundary condition (b.c.) needs to be supplied there. Given the 
smooth and, in particular, univalued nature of the eigen- functions (pujii^) there is a natural 
way to formulate a unique time evolution simply by following the usual recipe: decompose 
any incoming wave packet into (pi^ and then endow each component with an exp(iujt) time 
dependence. Since the eigen-functions (pujif) are smooth at r = 0, it is just natural in this 
formulation to consider both sides of the singularity to be connected and communicating. 

For negative-mass Schwarzschild the situation is different from the charged case, though 
our approach still provides a unique time evolution. In this case the r = singularity "flips" 
from being time-like at r < to space-like at r > 0, suggesting that gluing these two pieces 
would make little physical sense. Analyzing the wave equation suggests that indeed in 
the Schwarzschild case the field in the r < region evolves without any connection to the 
r > region. Unlike the RN case, here the eigen-functions (l)u){r) generically have a log 
singularity, hence they cannot be extended across the singularity in a univalued fashion. 
Yet, a single regular 01,^ (r) solution exists for each w. In this case a natural choice of 
boundary condition suggests itself - the so-called "regularity b.c." - i.e. the demand that 
(t>iuj{^) be bounded at the singularity. Since we have a single regular solution for each w, 
any wave packet coming from infinity will have a unique decomposition into these regular 
functions (/>ia;(r), and hence a unique time evolution. (This is to be contrasted with the RN 
case, in which there are two regular solutions for each w, and therefore, the time evolution 
depends on initial data from both sides of the singularity.) This "regularity b.c." can 
be physically interpreted as a reflecting b.c. While normally one can define reflecting b.c. 
either by Dirichlet or by Neumann b.c. (or a mix) the "regularity b.c." at Schwarzschild 's 
r = do not leave us such a choice ^ . 

Thus we find a natural way to define (scalar wave) physics in the presence of certain 
time-like singularities (by "defining physics" we mean a recipe that allows to predict the 
time evolution in spacetime) simply by passing from particle probes to waves (we call 
such singularities "wave-regular"). We discuss two mechanisms: for RN we continue the 
space-time across the singularity allowing for cross-communication, while for negative-mass 
Schwarzschild the "regularity b.c." are unique and define complete reflection. 

We would like to pause to make a few comments on this picture: 

• We consider this to be an important observation in the search for a resolution of black 
hole singularities, at least those of the time-like type. Yet, for this observation to 

''it was found already in by different methods that there exists a unique time evolution for negative- 
mass Schwarzschild, which is probably the same one we explicitly describe and physically interpret here. 
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be more than a mere curiosity it is required that several additional tests be satisfied 
such as regularity of higher-spin waves and regularity against various perturbations 
to the equations (some of which we perform). 

• At the classical level there are issues of interpretation, for instance, the existence of 
these spacetimes or others with the relevant singularities, the observers living in them 
and implications to Cosmic Censorship, which we discuss in section ^. 

• At the level of quantum gravity additional issues arise, including the validity of the 
wave equation and Hawking radiation whose discussion we defer to section |8| as well. 
Here we only want to stress that one should be cautious not to jump to conclusions 
and "legitimize" spacetimes with such singularities as possible solitons. 

After reviewing the RN and Schwarzschild backgrounds in section ||, we explain our 
main observation in section ^ We formulate the wave evolution for arbitrary wave-packets 
coming from infinity (i.e. from large negative r). However, the issue of specifying initial 
conditions on a spacelike hypersurface crossing the r = singularity becomes non-trivial 
and is relegated to section |^. 

Once the unique evolution of the field is formulated, we are in a position to investigate 
scattering phenomena. In section Q we study the physical properties of the RN singularity 
by computing the cross section for transmission. Even though the ODE for (jj^^ is regular, the 
effective potential diverges as T4ff ~ —c/r* ^ in the vicinity of the singularity, where r* = 
r*(r) is the canonical "tortoise" coordinate. Therefore, one would expect the singularity 
to have an effect even on incoming waves of very high energy - unlike a regular point in 
spacetime, which becomes "transparent" in this limit. Indeed we show that the high-energy 
transmission amplitude for fixed angular momentum I approaches 25% rather than 100%. 
In this respect the singularity behaves as a beam-splitter. Summing over / we find the high- 
energy dependence of the total cross section to be oc where uj is the frequency of the 
incoming plane wave. This behavior stands between the u:^ dependence of absorption for a 
black hole with well-defined area, and the dependence of scattering off an elementary 
particle in field theory. 

Generalizations of the backgrounds to higher dimensions and for the 2d black hole are 
presented in section ^. We determine the transmission amplitude and total cross-section as 
in the 4d case. The only significant difference is that while in 4d (and 2d) RN the variable- 
separated wave equation is regular, for c? > 4 it is singular-regular. Yet, all the solutions 
(j)uj{r) are still univalued for a continuation through the singularity. For Schwarzschild, on 
the other hand, there are no significant changes in passing from d = 4 (and d = 2) to d > 4. 
Hence, we also present the results for d = 4 Schwarzschild in section 5. 

Next we challenge our picture by adding various perturbations in section |6|. So far 
we discussed only probes which were minimally-coupled scalar fields. Here we preform 
several perturbations: first we consider a field with both mass and charge, then we add 
interactions, and finally we consider the effect of back-reaction. 

In section ^ we discuss some of the conditions for a general timelike singularity to be 
wave-regular, including the necessary constraint on initial conditions specified on a spacelike 
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hypersurface crossing the singularity, the uniqueness of decomposition of an incoming wave- 
packet, and the univalued property of the eigen-functions (pi^. In section ^ we summarize 
and discuss the issues mentioned above, and hst some open questions. Some technical 
details are given in the appendices. 

Note that, in principle, this study could have been carried out a long time ago and, 
indeed, it turns out that there were some attempts in this direction. However, it seems 
that none took the conceptual freedom to cross the singularity and that actually led several 
authors to pronounce RN to be wave-singular. Wald |^ gave a general discussion on 
extending the definition of the wave equation to include non globally-hyperbolic spacetimes 



relying on the mathematical notion of the "Priedrichs extension." Horowitz and Marolf |10| 
showed that some dilatonic black holes are wave-regular, where they consider a Schrodinger- 
like equation rather than the wave equation and demonstrate "essentially self-adjointness" 
of the associated operator. Ishibashi and Hosoya |^] were able to claim wave-regularity for 
negative-mass Schwarzschild by going back to the wave equation together with changing the 
metric in function space accordingly from the ordinary L^-metric to the Sobolev metric, but 
using the same criterion of "essentially self-adjointness." However, they do not comment on 
the physical nature of this resolution. Related work also appears in |11|. Finally, Jacobson 



|12U observed in passing the regularity of the time-separated wave equation near the RN 
singularity while studying the semiclassical decay of near-extremal black holes. 

Note added at second version: it is interesting to consider smooth approximations 
to the singular geometries. We found that if one smoothly approximates the RN singularity 
one is required to add a finite amount of negative energy matter, which is what one would 
expect due to the wormhole nature of the singularity region. The details of the calculation 
are given in section |2| in the part titled "Negative-energy content of the smeared singularity." 
The subsequent "perspectives on the r < region" part in section ^ as well as the discussion 
section should be read with this additional fact in mind. We gratefully acknowledge the 
valuable comments of G. Gibbons and J. Katz on this matter. See |l^ for a closely 
related case. 

2. Review of the Reissner-Nordstrom and Schwarzschild backgrounds 



Basic features of the Reissner-Nordstrom geometry 

The background we wish to consider is the Reissner-Nordstrom (RN) geometry 16 1. 
In Schwarzschild coordinates {t,r,9,ip) the metric is 

ds^ = -fdt^ + f-Ur^ + r^dn'^, (2.1) 

where 

and dn'^ is the two-sphere metric, 

dO^ = dO^ + sin^ 6 dip"^ . 
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The electromagnetic four-potential is given by ^ 



A = -^dt . 
r 

For Q = the metric reduces to the Schwarzschild metric [17|. 

For the (under-extremal) non-extremal black- hole case, M > \Q\ > / has two positive 
roots, located at 

r± = M± - Q2_ 

Both r = r_|_ and r = are null hyper surf aces, known as the event horizon and the inner 
horizon, respectively (in Schwarzschild there is only one such root, r = := 2M, which 
is the event horizon). Note that the roles of r and t as space- like and time- like coordinates 
interchange at each root of /. 

The line element ( ^.iD is singular at the horizons, but this is merely a coordinate 
singularity which may be removed by transforming to other coordinates (e.g. the Kruskal- 
Szekeres coordinates discussed below). On the other hand, r = is a true physical singu- 
larity (e.g. the Ricci scalar diverges in the RN case), and we focus on this singularity in 
the present paper. 

Usually one considers only the r > part of the RN geometry to be physical, but 
we will find it natural to consider the r < region as well. Note a crucial difference in 
this regard between the r = singularities in Reissner-Nordstrom and in Schwarzschild: in 
Reissner-Nordstrom / has a second-order pole at r = 0, 

/ = ^' (2-2) 

and the t coordinate is time-like on both sides of r = 0, while in Schwarzschild / has a 
first-order pole 

2M , , 

/ = — (2.3) 

and the coordinates "flip" their space- like/time- like nature: t is space-like for small and 
positive r and it is time- like for negative r. For this reason (as well as for other reasons), 
we will later find it natural to glue the positive-r and negative-r patches of Reissner- 
Nordstrom over the singularity, but not in the Schwarzschild case. The Penrose diagram of 
the maximally extended Reissner-Nordstrom spacetime, including the r < region, is given 
in figure ^. This is to be contrasted with figure ^ which provides the Penrose diagrams of 
the regions r > and r < in the Schwarzschild case. These diagrams - and particularly 
the fact that the singularity is spacelike at r > and timelike at r < - indicate the 
difficulty of matching the two regions at r = 0. 



Kruskal-Szekeres coordinates 

Although the Schwarzschild coordinates are perfectly suitable for describing the r < 
region and its matching to r > 0, insight into the "completed Penrose diagram" may be 

*The minus sign was entered in order to conform with the rest of our sign conventions and to produce 
( ^.19| ) correctly. 
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achieved by using the Kruskal-Szekeres coordinates, as we now describe. (In a first reading 
one may skip the text and only view the completed Penrose diagrams). 

The {r,t) coordinates (or "Schwarzschild coordinates") define several metric patches, 
which are separated on the r axis by r = 0, r_, r+ in Reissner- Nordstrom and r = 0, rg 
in Schwarzschild. In order to see how these patches, and especially the r < patch, are 



connected one transforms to the Kruskal-Szekeres coordinates [18, 19 1, and later performs 
a conformal transformation to obtain the Penrose diagram (in appendix ^ we review the 
procedure) . Here again we will find a qualitative difference between the Reissner-Nordstrom 
singularity and Schwarzschild. 

In Kruskal-Szekeres coordinates {U, V), one uses a function g(r) to implicitly define r 
through 

-UV = g{r). (2.4) 

This function is defined by 

g{r) := exp(±2Kr*) , (2.5) 
in terms of the "tortoise" coordinate r*: 

dr* := J . (2.6) 

K is the surface gravity of the relevant horizon and the choice of sign is correlated with the 
choice of sign for U and V. 

Let us look at the functions g{r). For RN g±{r) are given by (see figure |l]) 

«+ 

c/+(r) = ^-^ - 1^ ^-^-1^ exp(2K+r) 

g^r) = (l - (l - V^) exp(-2K_r) (2.7) 

where 5+ covers the "plus plane" which includes the outer horizon at r+ : r_ < r < 00, 
and g- covers the "minus plane" including the inner horizon at r_ : — 00 < r < r+ . From 
the form of g- one realizes that in these coordinates the r < region is naturally included 
and it lies at {—U- V-) > 1, which is adjacent to the < r < r_ or < —U- V- < 1 
region, and as such we include it in the Penrose diagram (figure ^). 
For Schwarzschild on the other hand, we have 

gscir«{r) = (^-^ - 1^ exp(r/ro). (2.8) 

One notes that unlike the g± which were monotonic gschwi''') has a minimum at 5'schw(^ = 
0) = —1 (see figure exactly because /schw('") changes sign at 0. Thus for —UV < — 1 
there is no solution for r while in the region —1<—UV<0 there are two solutions one 
of them with negative r. Hence, the r < region is naturally included and becomes a 
second cover over the inside of the black hole (0 < r < tq) which we incorporate into the 
Penrose diagram in figure ^. This second cover differs from Reissner-Nordstrom where the 
r < region was located side by side with the "ordinary" r > regions. This difference 
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Figure 1: The two functions g±{r) which determine r in Kruskal-Szekeres coordinates for Reissner- 
Nordstrom imphcitly through —U±V± ~ g±{r)- g+ is defined for the r > r_, namely on the 
plus plane, and g- is defined on r < r_|_ namely on the minus plane. Note that both functions 
are monotonic with range (—00, +00) and hence there is always a unique solution for r. For 
(— [/_ 1^-) > 1 on the minus plane r is negative. In drawing the figure the values r_ = .5, r+ = 1 
were used. 

is in tune with our later conclusion that for Reissner-Nordstrom one should glue the two 
regions over the r = singularity while in Schwarzschild they should be considered to be 
disconnected. 

Negative-energy content of the smeared singularity 

The electro- vacuum geometry constructed here has a curvature singularity at r = 0. 
One may consider modifying the metric function in the very neighborhood of r = 0, say 
at some |r| < ro, so as to make the geometry smooth. In this way one obtains a wormhole 
smoothly connecting the positive-r and the negative-r parts. However, any such smearing 
will imply the presence of additional energy- momentum in the region \r\ < rg, as would be 
determined by applying the Einstein operator to the modified metric. The specific form of 
the additional energy-momentum tensor will depend on the specific choice of smearing. It 
is well known, however, that any wormhole of this kind requires "exotic" matter fields, i.e. 
an energy-momentum source which violates the weak energy condition. A simple way to 
understand it is to consider a congruence of null rays propagating inward radially. In this 
setting Raychaudhuri's equation tells us that as long as the energy condition is satisfied 

log(yl) < 0, where A is the transverse area (of the two-sphere) and A is an affine 
parameter along the ray. Since A (and log(A)) is asymptotically decreasing in an inward 
radial motion, but after going through the "wormhole" it is again increasing, there must 
be an intermediate region where the weak energy condition is violated. 

We wish to evaluate the amount of negative energy required to support this wormhole- 
shaped smeared geometry, particularly at the limit ro 0. To be more specific, we consider 
here a thin-shell model ||2^. Namely, we directly match the geometry at r > ro to that 
at r < — ro, through a thin shell located at |r| = ro, for some < ro ^ Q. This thin 
shell is a timelike hypersurface parameterized by the three coordinates {t, 6, ip). Note that 
the induced three-geometry is continuous across this shell: g^e and g^pip only depend on r^, 
and gu (which originally has different values at the two sides) becomes continuous after a 
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Figure 2: The complete Penrose diagram for the (under-extremal \Q\ < \M\) Reissner-Nordstrom 
background including the regions with r < which are not drawn usually. These regions are glued 
to "the other side" of the time-like singularity, which is seen as naked from that side. The wavy 
lines denote a typical scattering experiment off the singularity. Waves from T~ are sent to the 
singularity. Some of them are reflected back to X+ while the remainder is transmitted to r > 0. 



trivial rescaling of the t coordinate at e.g. r < — ro. 



The shell's contribution to the energy- momentum tensor may be expressed as |21] 



T''''{x) = Sf'^x)6ir]) , (2.9) 

where S'^'^ is the shell's surface energy-momentum distribution, and r] is the proper-length 
parameter along radial lines of constant t, 9, ip, with r/ = at the shell and > at r > tq. 
The surface distribution S'^'^ is determined from the jump in the extrinsic curvature, 

S} = -i (W} - ^ . (2.10) 

Hereafter k run over the three coordinates (t, 9, ip), and 

{K)) = Kf+^ -Kf-\ (2.11) 
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Figure 3: The function .gschw('') which determines r in KruskahSzekeres coordinates (for the 
Schwarzschild black hole) implicitly through —UV — gschw{r)- Note that while for gschw > there 
is a unique solution for r, for — 1 < gschw < there are two solutions, one of them with negative r 
which lead to a double cover in Kruskal-Szekeres coordinates. Finally for gschw < — 1 there are no 
solutions. 




Figure 4: The complete Penrose diagram for a Schwarzschild black hole including the negative- 
mass regions which appear as a double cover of the inside of the black (or white) hole. The 
negative-mass regions are shown separately on the side, rotated by 90 degrees so that time flows in 
them upward as is conventional. Here we will find no communication through the singularity. 



where Kj is the shell's extrinsic curvature with respect to the geometries at r > rg and 
r < — ro, respectively. (We are using here General-Relativistic units, G = c = 1.) The 
extrinsic curvature may conveniently be expressed in terms of Israel's "natural coordinates" 
|2C], which we take here to be the three hypersurface coordinates t,6,Lp and the above 
proper-length coordinate ry: 



K. 



dg 



dr] 



(±) 



9, 



-1/2 



dr 



(±) 



(2.12) 



with the derivatives (and grr) evaluated at the (-I-) or (— ) sides of the shell, respectively. 

We are primarily interested here in the amount of (negative) energy contributed by 
the shell, represented by Sf: 



S 



■iTT 



{K% + {K)l 



Ait 



(2.13) 
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A straightforward application of Eq. ( 2.12| ) yields 



y^l-2M/r + Q2/r2' 



(2.14) 



r=±ro 



hence 

l^j^^e ^ yi - 2M/ro + QVr^ + ^1 + 2M/ro + Q^/r^ ^ ^^^^^^ 

Let p denote the shell's energy density as measured by a static observer: 

P = -T} = -j-{K)',5{v) (2.16) 

(recall that T/ is invariant to a rescaling of t). The total shell's energy is obtained by 
integrating over the 3-volume, i.e. over r] and over the shell's area: 

Eshell = 47rro j pdt] = -r^ {K)g 

= -ro (^^l-2M/ro + Q2/r2 + ^1 + 2M/ro + Q^/r^^ . 

Finally, we calculate Eg^eii at the limit where tq shrinks to zero: 

E'^j.^ii^Wui Eshell = -2Q. (2.17) 

ro^O 



Perspectives on the r < region 

The metric ( |2.lD is invariant under the transformation r — > — r, Af —M, and hence 
we shall often refer to the r < region as the "negative mass" region. To an observer in 
the asymptotically-flat region at r < the central object will appear as one of negative 
mass, namely, a gravitationally-repelling object.^ 

We alert the reader that negative masses have very unusual features: 

• A negative mass produces anti-gravity, namely it repels all masses. 

• As a result of Newton's second law F = ma, the acceleration of a negative mass is 
reversed to the force acting on it. 

• As an immediate result of the two properties above, if one were to place two bodies 
initially at rest, one with a negative mass and the other with a positive mass, both 
will accelerate in the same direction going from the negative mass to the positive 
one (furthermore, if the two masses are of the same magnitude, they will uniformly 
accelerate forever). 



^This is also the situation in the Kerr geometry: in the asymptotically-flat region r < the central 
object is gravitationally repelling. 
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Figure 5: The complete Penrose diagram for an over-extremal {\M\ < \Q\) Reissner-Nordstrom 
black hole. It is composed of a positive-mass RN glued over the singularity to a negative-mass RN. 
There are no horizons and the singularity is naked on both sides. 



Due to these unusual features (and possibly others) , negative masses are often excluded 
(however, see |2^). We shall not attempt here to provide a conclusive answer to this 
question. Instead, we will take a pragmatic approach of exploration by allowing negative- 
mass spaces in the current context and being alert to the appearance of a problem such as 
a resulting inconsistency. We will not find any such inconsistency in this paper, and so we 
regard this question as remaining open. 

In the (under extremal) RN black-hole case, M > \Q\ > 0, one can reach the r = 
singularity either from the negative-mass asymptotic region (i.e. from large negative r), or 
from r — > r_, the inner horizon. Looking at the Penrose diagram of the maximally extended 
spacetime, figure |2|, one may consider reaching the inner horizon by jumping into the black 
hole from the ordinary (positive-mass) asymptotic region. The inner horizon of the pure 
RN geometry is a perfectly-regular hypersurface. However, it is known already for some 



time that the inner horizon is unstable to perturbation and should become singular |23, 24| 
in "realistic" charged black holes. For this reason it was often believed that the inner 
horizon cannot be crossed. Recent investigations, however, indicated that the singularity 
at the inner horizon is weak, that is, the tidal forces are too weak to harm physical objects 



1 25, 26, 27 1 . Therefore, there is no obvious reason to exclude objects falling into the 
black hole and arriving at r = through the inner horizon. For conceptual simplicity, 
however, we shall primarily consider here observers arriving at r = from the negative-r 
asymptotically-flat region, namely, from beyond the singularity. 

In the case \Q\ > M (figure |5|) there are no horizons, but simply two asymptotically- 
flat regions, one with positive mass and the other with negative mass (in this case there is 
no reason to transform into Kruskal-Szekeres -like coordinates). In both regions there is a 
naked singularity at r = 0, and we glue the two patches there. 

Geodesies 

We end this section with a review of particle motion in this background, and a deriva- 
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tion of the effective potential for particles (or geodesies) which we later compare with the 
effective potential for waves. The geodesic equation for a particle of mass m can be derived 
from the Hamiltonian system 

n = ^g'"'p^p, = ~m'' , (2.18) 

where the Hamiltonian is automatically with respect to a multiple of the proper time r. 
For a particle of charge q one replaces p ^ p — q A. 

Since H is independent of t and cf) we have the constants of the motion E := —pt, I := 
P(f,. The motion is then described by 

rv • 2 I T ^geodesic / \ 

= r + (r) , 

^^g^eodesic(^) = + m^) fir) -(^E-^-^y . (2.19) 

In particular, around r ~ 0, VcS ~ for Z > [and V^cff ("^-^ — Q^)^ foi" ^ = 0]. Note 
that the effective potential given here does not have the ordinary units of potential, but it 
has the advantage of being general and valid both for massive and massless particles. In 
order to get the physical potential in a Newtonian approximation one needs to rescale V^ff 
by an appropriate power of m. 

3. Waves Eire smooth 
3.1 Set-up 

A spacetime is usually defined to be singular when it is not geodesically complete, namely 
with respect to point-particle probes. Yet, all the fundamental interactions are formulated 
for field theories, where particles are derived concepts constructed from appropriate wave 
packets. Moreover, the mashing of Quantum Mechanics and Special Relativity requires the 
quantum theory of fields. Thus it is natural to probe spacetimes with waves rather than 
point-particles. 

Here we will confine ourselves, as a first step, to scalar fields. Initially we will further 
assume a minimally-coupled free field with no mass or charge, and later these assumptions 
will be relaxed and shown not to affect the results. 

When should we consider a spacetime to be wave-regular? A central objective of 
physics is to provide predictions for the results of experiments, or a well-defined time 
evolution, and it is exactly this predictability which is at risk in the presence of a singularity. 
We would therefore consider a spacetime to be wave-regular if the time evolution is unique 
for any physically acceptable initial conditions (to be discussed below). Technically the 
issue is whether one can supply boundary conditions at the singularity such that the time 
evolution is regular and whether these are unique. 

In "normal" (i.e. non-singular) situations there are two different set-ups for initial 
conditions: (i) the Cauchy formulation, in which the field and its time derivative are 
specified on a spacelike hypersurface; and (ii) the characteristic formulation, in which 
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the field is specified on two null hypersurfaces. This second formulation has a special, 
widely used variant (iia) in which the null hypersurfaces are taken to be at the null past 
boundaries of the region of spacetime under consideration - e.g. past null infinity (for an 
asymptotically- flat region) , and a past horizon (for an asymptotically-flat region outside a 
black hole). 

In our problem - formulating the time evolution of the field at the two sides of the r = 
singularity - the Cauchy formulation (i) has a potential problem, because any spacelike 
initial hypersurface must cross the r = singularity. It is a priori unclear what are 
exactly the "physically-acceptable initial conditions" near r = 0. On the other hand, the 
characteristic formulation (iia) is free of this potential problem. The two null hypersurfaces 
which naturally suggest themselves for this formulation are past null infinity (X^) of the 
negative-r universe, and the inner horizon r_. These two hypersurfaces nowhere intersect 
the singularity, or even get close to it. 

We shall therefore base our analysis of time evolution on this characteristic formulation. 
That is, we probe the spacetime by wave-packets prepared and sent in from the asymptotic 
regions far away from the singularity. Indeed, it is physically reasonable to characterize 
the singularity by its response to any possible scattering experiment. Motivated by these 
considerations, we shall define a spacetime (or a region of spacetime) to be wave-regular if 
the time evolution is unique for generic initial conditions specified at the relevant past null 
boundaries of this (region of) spacetime. 

While the characteristic formulation (ii) avoids the issue of putting constraints on initial 
conditions near the singularity, which is indeed non-trivial in hindsight, it does require a 
global point of view, whereas our resolution is essentially local. This local property is 
better understood in the Cauchy formulation (i) and hence we will briefiy review both 
formulations in this section. Actually one could define also a hybrid of these methods 
which keeps both locality and some isolation from the singularity: namely, a formulation 
(ia) in which the field and its derivatives are supplied on a spacelike hypersurface but are 
required to vanish in a neighborhood of the singularity. ^ 

In what follows we shall show that both the — oo < r < r_ piece of the RN spacetime 
and negative-mass Schwarzschild are indeed wave-regular. We shall define first the time 
evolution in a straightforward manner in the Cauchy formulation, and then after defining 
a canonical form for the radial equation and the effective potential we shall state the char- 
acteristic formulation. A treatment of the interesting subtleties in the Cauchy formulation 
will be deferred to section 0. 

3.2 Separation of variables 

Let us proceed by writing down the wave equation for a minimally coupled massless scalar 
field 4> in the RN (or Schwarzschild) background ( p.ll ). After separating the angular vari- 



This variant is slightly more restrictive than just a re-formulation of the initial- value setup, as it demands 
that the field strictly vanishes at some open neighborhood. Nevertheless it is still sufficiently general to 
allow for physically-meaningful scattering experiments. 



ables 





- 14 - 



where Q, denotes the angular variables ^ {6, (p), one has 

= 00= \-r'du + \dr.fr^dr-^-^^ 



(3.2) 



The wave equation (|3.2| ) is singular at r = 0: after multiplying the equation by / to extract 
dtti the part containing e.g. the second-order r derivative becomes f'^ d„i and diverges 
as r~^. This is not surprising since we know that r = is a curvature singularity. 
Still, we may attempt a straight-forward separation of the time variable, 



getting 



= 



4>iir,t) = (l)^i{r)exp{iuJt) , 



(3.3) 



(3.4) 



It will be convenient to normalize O^ji in two different ways: if we want to have a unit 
coefficient for drr then we work with O^^i while / D^i achieves a unit coefficient for a;^. 
Indeed, writing down D^i explicitly, we have 







/ ^ OojI 4>ojI = 



(3.5) 



and at this point we need to discuss the Reissner-Nordstrom case and the Schwarzschild 
case separately. 



Reissner-Nordstrom 

Since Jrn = is regular at r = 0, the variable-separated equation ( |3.5[ ) is surpris- 
ingly completely regular. This observation is at the heart of our proposed resolution. To 
leading order it is given by 







(1 + . . .) drr + 0(r°) dr 



;i + ...)^ + (i + ...)^-^ 



(3.6) 



where ellipses always denote corrections which are higher order in r, and from now we will 
often write (j)^^ instead of cp^^i suppressing the index I for clarity. Note that the equation 
is not only non-singular but also analytic (namely, 0aj,T-r is an analytic function of (p^j, 4'uj,r 
and r) at r = 0, for any I and io. 

Since all the solutions for this equation are smooth at r = 0, a natural way to define 
the time evolution suggests itself: given an initial wave packet (which propagates from 
either 2~ and/or the inner horizon) decompose it into radial eigen- functions (j)uj{r) . Since 
the eigen- functions (pujif) are smooth at r = 0, the time evolution defined above will be 



^The notation fl for the angular variables is convenient for the generalization to higher dimensions where 
there are more angular variables and there is no need to explicitly define them. 
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smooth there. Consequently, parts of the wave packet will cross this point. It is therefore 
essential to include both sides of the r = singularity in our analysis. 

In equations: given initial conditions (/'i(r) := (j){r,t = 0) and (pi^r) := i;^>(r, i = 0) (we 
chose here the initial hypersurface t = without loss of generality), one determines the 
decomposition (j)^{k), (j)-{k) that satisfies 



Mr) = J {Mk)+^-{k)) Mr)dk 

Mr) = [iuj iMk)-(p-ik)) Mr)dk (3.7) 



where here uj := \k\, and c/i^ denotes two independent solutions (one for positive k and one 
for negative k) of the wave equation O^i = (with a; = |fc|). * Now the time evolution is 
defined naturally to be ^ 



0(r,t)= J {((>^{k)e-xp{+iujt) + (f)_{k)ex-p{-iujt))((>k{r)dk . (3.8) 

In the characteristic formulation the decomposition into modes may be done by Fourier- 
decomposing the waves coming from I~ or from the inner horizon in exp(i uj v) or exp(z a; u), 
respectively, where u and v are the two Eddington-like null coordinates, and associating a 
certain radial function (f)uj{r) with each Fourier component. This will be described in more 
detail in subsection |3.4 

At first it may look as if the singular partial differential equation has turned 

completely regular. Later we will see, however, that some imprint of the singularity remains. 
This is expressed both in the unusual large-w limit of the transmission amplitude through 
r = (see section ^) and in the constraints for Cauchy initial conditions at r ~ (see 
section 0) . 

Schwarzschild 



Using the leading behavior of /schw at r = (|2.3D , the variable-separated equation 
|3^ ) is given to leading order by 







(3.9) 



One notices that this equation is not regular anymore, but rather the dr and the /(/ + 1) 
terms are singular. Actually it falls into the so called "regular-singular" category of mild 
singularities (see appendix^ for a rudimentary review). 

Looking at Eq. ( |3.9| ) one realizes that the characteristic exponents (defined by the 
attempted leading-order solution cp^ = r^) are pi = P2 = 0, and thus the two independent 
solutions behave as cpi^i ~ r^, 02w ~ log(r) at r ~ 0. So unlike Reissner-Nordstrom one 
of the solutions is singular, and moreover does not have a single- valued continuation to 



4>k are usually chosen such that (jik — exp(i kr* )/r at large r* (large negative r). 
^This amounts to taking the radial function 4>^i{r) of Eq. (B.3) to be 4>+{k = + 4>+{k — —u)) for 
UJ > Q and <j)-{k — +ui) + (j}-{k — —uj) for cj < 0. 
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negative r, and there is no sense in gluing together positive r to negative r. It is perhaps 
fortunate that this is the case, since gluing the two sides would require dealing with the 
"flip" in the nature of (r, t) from space-like to time-like (and an arbitrary choice of the 
direction of the arrow of time). Still we find it remarkable that one of the solutions does 
turn out to be regular. This half-regularity allows a natural definition of the time evolution, 
as we shortly describe, where this time there is no transmission through r = 0. 

The situation here is analogous to the relation between a free field (j) on the infinite 
line — oo < r < +00 to a free field on the semi-infinite line < r < +00, which could be 
e.g. the radial coordinate in some higher dimension, with Dirichlet or Neumann boundary 
condition at r = 0. For the infinite line there are two </>(^ solutions for each uj, which allow 
a decomposition of wave packets arriving simultaneously from both —00 and +00. On the 
other hand, for the semi-infinite line the boundary condition selects a single <j)i_j solution for 
each CO. This still allows for a decomposition of a wave packet, which this time can arrive 
only from +00. The situation for Reissner-Nordstrom is similar to the infinite line, while 
negative-mass Schwarzschild is similar to the semi-infinite line (we consider the r < side 
of the Schwarzschild singularity since there t is time-like and one can consider scattering 
experiments). 

Thus we naturally define the time evolution for negative-mass Schwarzschild by de- 
composing the initial conditions (defined only at r = —00) into the regular eigen-functions 
(pii^ alone and then proceeding as usual to define the time evolution. Since (pii^j is smooth at 
r = 0, the time evolution will be smooth there too. As it turns out it was found already in 
Q by somewhat formal arguments that there exists a unique time evolution for negative- 
mass Schwarzschild, which is probably the same as the one we just explicitly described. 
For RN on the other hand we differ as we get wave-regularity after considering both sides 
of the singularity, while they do not. 

Perspectives on the boundary conditions for Schwarzschild 

In the Schwarzschild case we simply impose the boundary condition that the scalar 
field be regular at r = 0. This does not conform with either the Dirichlet or Neumann 
b.c; namely, neither (f) nor its radial derivative vanishes at r = 0. In fact, the Dirichlet 
or Neumann b.c. are both inconsistent with the asymptotic behavior near the singularity 
(due to the presence of a divergent mode). 

Yet, there is a physically appealing procedure to obtain our regularity b.c. from the 
Dirichlet or Neumann b.c. (or a mixture). We can cutoff space at some ri < close to the 
singularity, and then take the limit ri 0. Since all modes are well-behaved at r = ri, 
we are free to impose either Dirichlet or Neumann b.c. there. Let us denote as before 
the solution regular at r = by and some other independent and necessarily singular 
solution by <j)2Lu- We are looking for a linear combination of the two, </>^ = acpi^ + b(j)2uj, 
that will satisfy the boundary conditions at ri. If we impose Dirichlet we need 

= Mn) = a(l)iUri) + 602c.(ri) . (3.10) 

As we take ri — > c^i^j remains finite while 4)2^ diverges with a log singularity, and hence 
we must take 6 = 0, namely the "regular solution b.c." 
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If we choose Neumann conditions instead, we note that again while remains finite 
(j)'2^^ diverges with a 1/r singularity and hence we are again forced to take 5 = 0. Altogether 
we find that in both cases as ri — > we are left with a unique b.c, namely our "regular 
b.c." 

3.3 Normal form of the radial equation 

Transforming to a new radial coordinate (the so called "tortoise coordinate") r*(r) defined 

by 

dr* = (3.11) 

fir) 

and to a new field variable 

(j)* := (3.12) 



(and correspondingly (/>* := r(f)i_j), the wave equation (|3.2|) and the radial equation ( p^ ) 
take the standard forms 



r' [dtt - d„r. + V^sir*)] <P* = (3.13) 



and 



[-d„r* + Vesir*) - u;2] = 0. (3.14) 

Here the effective potential Ves is given by the effective potential for particles, 

with the substitutions P — > /(/ + 1) and E ^ 0, plus an extra term AVcs'- 

V,s{r*) = AVes + / (3.15) 

where 

AV,s{r*) = - dr f . (3.16) 

r 

We find it instructive to explore the leading terms in the radial equation at r = 0, 
when this equation is expressed in terms of r* rather than r. 

We start by looking at the Reissner-Nordstrom case. To leading order near r ~ 



and hence 



(3.17) 



n4 2 

Ves = -2-^ ^ ^ . (3.18) 

Thus even though the wave equation in the form is smooth, the effective potential 
in the normal form is unbounded from below; hence we will study later (in section ^) its 
response to high-frequency waves, where all features of the potential other than its leading 
term oc r* ~^ are irrelevant. Note that such potentials of the type V = —c/r* ^ are special 
(see for instance [^). Considering a larger class of potentials V = —c/r* " with c, a > 0, 



^°in principle E is to be now replaced by uj, but the latter already appears explicitly in Eq. (3.14). Recall 
also that we assume at this stage a neutral massless scalar field, hence m = g = 0. 



-18- 



then a = 2 is a critical value describing a conformal quantum mechanics namely c 
is dimensionless. For < a < 2 the spectrum is bounded from below even though the 
classical potential is not, while for q > 2 the spectrum is not bounded from below. For the 
critical value a = 2 there is a critical value for c, c = 1/4. For c < 1/4 the potential is still 
mild enough to have a spectrum bounded from below (actually at zero), while for c > 1/4 
it is not. So the numerical coefficient c = 2/9 in ( |3.18| ) is important. We will see later (in 
section |5|) that c varies over our examples (it depends on the number of dimensions, but 
always c < 1/4 in the charged case), and that it determines the high-energy transmission 
properties of the singularity. 

At r* = the radial equation ( p. 14 ) has a regular singularity, and we may look for the 



characteristic exponents p defined by (p"^ = r* p. These satisfy the characteristic equation 



leading to the two roots 



p{p-l) + l = 0, 



12 , , 

Pi,2 = 3, 3, (3.19) 



which are independent of uj. Translating back from r* to r using the transformation 
( 3.12 , 07|) we get the exponents r'^and r^, exactly as expected for a smooth second-order 



equation. 

For Schwarzschild one may repeat the calculations based on the leading behavior of 



/schw (|]3|) and find 

„2 



r* ^ (3.20) 
4M' ^ ^ 



leading to 



eff = -4— = -T-77- (3-21) 



Thus we have here the critical value c = 1/4 (independent of the spacetime dimension, see 
section ^) . 

Attempting (/>* = r* ^ one obtains the characteristic equation 



which yields 



p{p-l) + \ = 0, 



Pi = P2 = -. (3.22) 



The equality of the two exponents implies the presence of a log term in one of these 
solutions. Again translating back the characteristic exponents to the r coordinate we get 
the r^, log(r) behavior which was derived above from Eq. 



if we replace cu co in (3.14) 
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3.4 Characteristic formulation 



Having defined the "tortoise" coordinate r* ( p.ll ) and the effective potential ( |3.15| ) we 



can now return to spell out in detail the definition of time evolution in the characteristic 
formulation. 

The coordinate r* approaches — oo as r — > -co and +00 as r — > r_. Notice that the 
effective potential decays as r"^ ~ r* at large negative r, and as / oc r — r_ (which 
is exponentially small in r*) at the inner horizon. Therefore, in each of these asymptotic 
boundaries the radial function is a superposition of the two asymptotic solutions 

^ exp(+za;r*) and ^ exp(-za;r*) . (3.23) 

Defining the Eddington-like coordinates 

V := t — r* , u:=t + r*, 

the above two asymptotic solutions at both edges r* ±00 read, for a particular u 
component of (j)* , 

(j)* = exp(icjn) and (p* = exp(iujv) , 
respectively. The non-decomposed field at both asymptotic edges takes the form 

where f±{v) and g±{u) are arbitrary functions (the "it" stands for the two asymptotic 
regions r* — > ±00, i.e. to the inner horizon and to large negative r). In the characteristic 
initial- value formulation one needs to specify the ingoing component f-{v) at X~ and the 
component g+{u) coming from the inner horizon. These two functions should uniquely 
determine the time evolution in the relevant piece of spacetime. This evolution is defined 
as follows. First, Fourier-decompose these two initial functions: 



f_[v) = j f^{u;v) exp{iu; yv)du;v, 
g+{u) = I gu{uJu) ex.p{iuj uu)duJu. 



Next, we define a set of basis solutions (j)*^^ Q^(r*,t) and '/'(o (^°^ each ujy and Uu) 

as follows: is the solution evolving from the characteristic initial conditions 

f-{v) =e-x.p{iujyv), g+{u)=0. 

Similarly, (p'^Q ^ ^ is the solution evolving from the characteristic initial conditions 

/_(v) = 0, g+{u) = ex.p{iu; uu). 

The two functions (j)*^^ and (p'^Q ^ ^ are constructed from the radial functions 4>lj{r*) (with 
LO = LOy and io = lOu, respectively); see appendix |C| for details. 
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Once (j)*^ and ^ are defined, the time evolution of the field (p* emerging from 
the characteristic initial data f-{v) and g+{u) is simply given by 

We conclude that the unique extension of the radial functions 4'lj{r*) across the r = 
singularity naturally leads to a unique time evolution for any set of characteristic initial 
functions f-{v) and g+{u). 



4. Transmission cross section in Reissner-Nordstrom 

Having demonstrated a resolution of the r = singularity we turn in this section to 
examine its nature through scattering. Here we shall only consider the charged (Reissner- 
Nordstrom) case. Full reflection from the negative-mass Schwarzschild singularity will be 



considered together with the d > 4 cases at the end of subsection |5.1| . We focus here on 
the transmission of high-frequency incoming waves through the r = singularity (from 
now on we will use "high-energy" to mean "high-frequency" by a slight abuse of language 
originating from the quantum theory). 

We first study the effective potential Ves, and especially its asymptotic form at small r. 
Then we compute the transmission and reflection coefficients for fixed / in the high-energy 
limit. Finally we combine that result with the high-/ limit, where the partial cross section 
vanishes. We sum over I and obtain the uj dependence of the (high-energy) total cross 
section for transmission. 

4.1 General features of the effective potential 

Consider first the qualitative behavior of the effective potential V^g. For concreteness we 
shall consider here a wave propagating from past null infinity of the negative-r universe 
towards r = 0. At large |r|, we have the usual centrifugal barrier 



^eff 



/(/ + !) 



(|r|>M). (4.1) 



At small |r|, the potential takes the asymptotic form 



,,s = Kl + l)^-2^ (|r|«Q). (4.2) 



Here we only keep the leading order in r/Q. We do keep, however, the term oc in the 
right-hand side because it is proportional to /(/ + 1). This term will become important 
for partial waves of sufficiently large I, which are involved in the calculation of the total 
high-energy transmission cross-section (see below). 

The contribution of any particular I to the total high-energy cross-section is suppressed 
by the kinematic factor k~'^ in Eq. ( 4.14| ) below (while it will turn out that a ^ k~'^). 
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Therefore, at high energy the total cross section is dominated by partial waves with Z ^ 1, 
which we shall assume hereafter. Such waves of relatively small or moderate energy will be 
reflected already at the large-|r| centrifugal barrier ( [4.1[ ). Only waves with sufficiently high 
energy, to > l/Q, will penetrate into the central region \r\ <^ Q. The scattering features of 
these high-energy waves will be dominated by the small-|r| potential ( |4.2| ). This potential 
has a peak value (hereafter we often replace l{l + 1) by which is justified because / 1) 



V — 

''max 



/6 



27Q2' 
located at 

\r\ = Tpeak ■= VSQ/l . 

It acts as a repulsive barrier (oc r~^) at Q S> I?"! > Vp^ak, and as a potential well (oc —r~^) 
at \r\ < Tp^n^. 

The above discussion immediately suggests that partial waves with lo <C ujpeak: where 

^3 

'^peak ■= \/Vtnax = /— „ , (4.3) 

will be fully reflected by the potential barrier at \r\ > Vp^ak- However, partial waves with 
uj » ojpeak will predominantly feel the potential well oc — at |r| < rp^ak- These waves 
will be partially transmitted, as we analyze in the next subsection. 

For later convenience we transform the small-|r| potential from r to r*, using = 

^ QlT^^ - ^ 

The peak is located at 

\r*\=r;,ak ■.= V3Q/l'. 
Although we are primarily considering here a wave coming from the negative-r past null 
infinity, the case of a wave propagating from r_ towards r = may be treated in exactly the 
same manner. Note that the small-|r| potential (4.2) or (4.4) is even in r or r*, respectively. 



Consequently the high-energy transmission amplitude for partial waves is the same for both 
cases. (Obviously the transmission probability is exactly the same, even for finite w.) 

4.2 Amplitudes at high energy and fixed I 

It is instructive to look at the scattering in the high-energy limit. Normally, the potential 
is bounded, hence the transmission amplitude approaches 1 and the reflection vanishes at 
this limit. However, here we have an unbounded (attractive) potential near r = r* = 0, 
which leads to partial reflection even at the high-energy limit. 

For fixed I and high energy, uo ^ ujpeak, the potential is dominated by 

^eff = =Kfr,sing- 



We therefore need to solve the scattering problem 



= [-dr^r* + Kff,smg(r*) - Uj'^W . (4.5) 
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Passing to the dimensionless variable 

X = —r* to, (4-6) 



the equation is transformed into 



Its general solution is 



0=[-d,,--^-l]r. (4.7) 



b* =Ci^l+ C2 Ct>*2 , (4. 



with 

(At = VxJi/q{x), 

^2 = J-l/6{^) (4.9) 

(for x > 0), where J„ are Bessel functions of order n (see appendix^). The coefficients 
ci , C2 will be determined below. 

In order to extend the solution (|4.8| ) across x = 0, we first express it in terms of r. 
Recall that at small r, 

^* ~ ^ 



3Q2' 

and, furthermore, r*(r) is analytic at r = 0. From the local behavior of the Bessel functions 
J^(A) near A = (see appendix Q), and that of r* near r = 0, it follows that <j)* and (j)2 
take the form 

<t>l = rAi,^,{r'), (4.10) 

where ^^-^yg(A) are some real functions (for real A), analytic at A = 0. Obviously both 
(/)t(r) and (t)2{f) have a unique real analytic extension through r = 0.^^ The form of Eq. 
( 4.10[ ) ensures that c/ij extends as an even function, and (1)2 as an odd one. Hence for 



negative x we get 

<l^i{x) = y/\x\Ji/&i\A)^ 

ct>*2{x) = -V\AJ-i/M)- (4.11) 

In order to get the coefficients ci , C2 - and the transmission and reflection coefficients 
- we impose the large-r* boundary conditions. Since we are considering here a wave 
propagating from (say) past null infinity (large negative r*), we assume that no waves are 
entering from the inner horizon (r* +oo). This corresponds to pure asymptotic behavior 



^^The negative sign in the definition of x was introduced in order for the incoming wave to propagate 
from large positive x values (null infinity) towards smaller x values. 
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Recall that 0i(r) and (j)2(r) must extend to r > as real and analytic functions, because the radial 



equation (3.5) is real and analytic (and 'j)i,ff)2 are real at r < 0). 
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(4.12) 



oc exp(— r*) = exp(i x) at x ^ — oo. Thus, in terms of x the asymptotic behavior at the 
boundaries is 

4>* = exp(ix) + R{u)) exp(— ix) at a; ^ +00 , 
(p* = T(u;)exp(ix) at x — > —00 . 

Note that since u; got scaled out of the equation T{uj), R{uj) are actually independent of 

UJ. 

The asymptotic form of the Bessel functions for large (positive) argument is given in 
Eq. (p^). For n = 1/6 it reads 



I) ^ V7r|x|/2exp(±i7r/3) (1 ± iV3)Ji/e{\x\) =F 2iJ_ 



-1/6 ( 



Applying it first to the inner-horizon boundary (x 
finds 



-00 



I, along with Eq. ( 4.11| ), one 



ci = ^/^^/2exp{-i^T/3){l - iV3)T{uj) 
C2 = —iV27rexp{—iTT/3)T{uj). 



Similarly, when applied to the null- infinity boundary (x — > -|-oo), along with Eq. (^), it 
yields 



cx = ^JVj2 exp(i7r/3)(l -hi\/3) + exp(-i7r/3)(l -ii\/3)ii(u;) , 
C2 = — i\/27r [exp(i7r/3) — exp(— z7r/3)i?(Lij)] . 

Combining the two expressions for c\ , C2 and solving for R and T, one obtains 

i? = ^i, T=\. (4.13) 

We conclude that for any /, at the high-energy limit a fraction \R}^ = 3/4 of the influx is 
reflected and a fraction |Tp = 1/4 is transmitted to the other side of the r = singularity. 
As was mentioned above, this applies to both waves coming from past null inflnity (of 
the negative-r universe beyond the singularity) and waves coming from the inner horizon 
towards the r = singularity. 

4.3 Summing over I 

The total cross section for transmission is given by a sum over all partial waves /: 

00 

z=o 

where k is the momentum of the incoming plane wave {k = lo for a massless held, which 
we assume from now on in this subsection). We shall consider here the cross section at the 
high-energy limit. As mentioned above, at this limit the total cross section is dominated 
by the contribution from modes / » 1, because each particular I is suppressed by the 
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factor uj ^, while it will turn out that a ^ uj ^. The result of the previous subsection, 
|T/|2 ^ 1/4, holds for I values that are not too large. However, for any fixed and large 
u) there is some large enough I where the transmission practically vanishes. For instance, 
for classical flat-space scattering off a rigid target of length scale L ("rigid" means here 
w-independent L), if the impact parameter b is much larger than L then the absorption 
vanishes. Since b = l/p (and p = uj in our massless case) we see that for / > Liv the 
absorption vanishes. We will now show that in our case the suppression of transmission 
already happens for / > {Qlo y^'^. 

A rough estimate of the total large-w cross section can be easily done based on the 
qualitative features of Ves at small r. For uj » ojpeak, where ojpeak is defined in (^), 
the above high-energy result |r,p ^ 1/4 holds. For u! significantly smaller than uip^ak 

the 

incoming wave will encounter the potential barrier 

^efr = g2/3(3,.)4/3 (4-15) 

and will be tunnelling suppressed there. [In fact, an effective tunnelling suppression requires 
that is greater than 1, where A is the tunnelling length scale in the barrier in terms 
of r*, i.e. from ( fl.lS ) A ~ {uj'^Q /l'^)^^/'^ . But this is automatically satisfied because 



wA ~ {u)Q is roughly {ujpeak/'^Y^'^ ^ which is S> 1 in the present discussion.] Thus, 
for a rough estimate of the large-w cross section we may substitute \Ti\^ = 1/4 for / <C Ic 
and negligible jT^p for I ^ 1^ where 

■= VsiQuj)^/^ (4.16) 



can be read off the expression for ujpeak ( |4.3| ). 

Equation ( ^.14| ) then yields the total cross section: 

a-h-^. (4.17) 

The total cross section may also be obtained from Ic by a simple argument: the critical 
maximal impact parameter for transmission is 6max = ^c/^, and hence the total size (cross- 
section) of the target for transmission is 

"max ^2 ^4/3- V^-^'^) 

In the above estimate we have not taken into account the contribution coming from 



the range / ~ /c- A more rigorous derivation of the cross section ( 4.17 , 4.18 ), which properly 
treats this range as well, may be done by analyzing the rescaling properties of the small-r 
scattering problem, i.e. 

{-dr- r* + Kff - = (4.19) 

with the potential ( |4.4| ). Transforming to a new dimensionless variable 

X := --^r*, (4.20) 
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the equation becomes 

[-d^^ + V{x)-n'']<p = Q, 

where 



9x2 (3^)4/3 

The boundary conditions at large |x| are, in accordance with Eq. ( [4.12| ), 

(p* ^ exp(ir2x) + R{Q) exp{—iflx) at x — > +oo , ,^ ^-j^x 

Hence in this limit of large to and large I the transmission amplitude depends only on 17, 

T = t(^). (4.22) 

The cross section is 

where we used / ^ 1 to replace 21 + 1 hy 21 and the sum by an integral. Finally, trans- 
forming from / to / := 1/{ijjQY/^ we find 

a-c^, (4.23) 



where 



c = 2^ j i\T{i-^)\^dl 



in agreement with the above estimate (4.17) 



Having found the total cross section for transmission, one may inquire about the total 
cross section for scattering. At large distances this is just like a Coulomb problem with a 
length scale M. Due to the infinite range of the gravitational force, we expect the cross 
section for scattering to be infinite for a massive field. However, for a massless field the 
scattering potential decays (after subtracting the standard centrifugal piece l{l + l)/r* 2) 
as ln(r*/Af)/r* 3. This is a short-range potential, hence we expect a finite scattering 
cross-section in the massless case. 

The peculiar uj dependence of the transmission cross section ( [4.17| ) through the singu- 
larity can be compared with two other systems. For black holes with non-zero area the large 
UJ absorption cross section through the event horizon tends to the geodesic cross section, 
and hence a ~ oj^ . For elementary particles, on the other hand, the high-energy scattering 
cross section cxscat ~ "^"^ which is the kinematical factor which multiplies the dimensionless 
scattering amplitude. So we see that this system lies between the two examples above. 

Another interpretation issue is whether the partial transmission/partial scattering be- 
havior at large iv should be interpreted as a beam splitter for single particles. The question 
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is whether a wave packet can be made small enough so that it fits into a target of size 
~ u;"^/^ - otherwise it will not be fully split. Imagine that we set a "lens" (one which 
affects (f) waves) at a distance f ^ uj^ ^ Q from the black hole. We know that the spot 
size at focus can be made as small as Xf/D, where A = 2tt/uj is the wavelength and D 
is the aperture size, and we here imagine that space is flat and we need to hit a target 
with the same cross section as the r = singularity. From Xf/D <^ uj~'^^^ we find that 
we need D > 

So such focusing can be readily achieved, even with a small size lens. 
This suggest that the singularity may indeed act as a beam splitter. Note, however, that 
in the above analysis we only considered incoming plane waves (or pure partial waves with 
well-defined I values). In order to verify the beam-splitter behavior one needs to explicitly 
analyze the transmission amplitude in the geometrical configuration considered above (i.e. 
a narrow beem focused on a small region near r = 0). This calculation is beyond the scope 
of the present paper. 



5. Generalizing the backgrounds 
5.1 Higher dimensions 

The background RN metric and gauge field in d > 4 dimensions are |29] 

A 



Q 



r-d—3 



dt 



The mass and charge are 



/ 



M 



r-d—Z 



y.2 (d-3) 



(5.1) 



{d-2) ^d-2^^d^-3 ^^d_-3^ 



16 vr 



Q = const y/r^ r_ 



d-3 



(5.2) 



where il.d-2 is the area of the d — 2 sphere, and we do not fix the constant in front of the 
charge which depends on the normalization chosen for the vector field. In order to get 
d-dimensional Schwarzschild one sets r_ =0, = tq. 

After separation of the time variable {du — > — w^) the scalar wave equation becomes 



1 



fr 



d-2 



drfr'^-^d, 



l{l + d-3) 

77^ 



. 



(5.3) 



Considering the charged case first, and comparing with the 4d case ( p.6| ) one notes a 
change: the equation ( |5.3| ) is not regular anymore at r = (since f r'^^'^ is has a pole). 
Actually it is regular-singular with leading behavior 



dr + l{l + d-2,)0 fr2(°'-4)) + u;2 Q 



,4(d-3) 



= . 
(5.4) 
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The two solutions of the characteristic equation are r", r'^^^. Since the difference of the 
exponents is integral it is possible a priori that the r^ solution has a log piece of the 
following form (j)^^ = r^ + . . . + log(r) r'^~^ + . . ., however it is seen not to be the case by 
explicitly Taylor expanding cp^ in the differential equation ( |5.4| ). 

For Schwarzschild, on the other hand, one finds the same leading behavior as in 4d 



drr + -dr + 

r 



(5.5) 



and thus the leading behavior of the two independent solutions is still r^, log(r). 

As in subsection we may find the normal form of the radial equation ( ^.31) . The r* 
coordinate is defined by 

dr* = j , (5.6) 

the field is redefined by 

0* := Ad-2)/2^ ^ (5^7) 
and the resulting effective potential is 

l{l + d-3) 



Aye. = ^^9.(r(-^)/V) (5.8) 



where we use the same notations as in subsection 3.3 



We start by looking at the charged case. The leading behavior near the singularity (at 
r = 0) is 



(2d -5) (r+ r_)('^-3) 
{d - 2) (3d - 8) 1 



4 (2d -5)2 r 



* 2 



(5.9) 



Given an effective potential Ves ~ —c/r* ^ the order n of the Bessel functions which appear 
in the high uj solution is given by (see appendix ^) 

= J - c . (5.10) 



Hence, the second equation in (|5.9|) determines 

, , 1 d-3 

Hkn = (5.11) 

generalizing the 4d result |n| = 1/6 ([4.91) . 

We note in passing that for the 2d black hole we will show in the next subsection that 
Vefi — — 3/(16r* 2) and hence one can define an "effective" or "equivalent" RN dimension 
to be dcs = oo or 11/4, and \n\2d = 1/4. 
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The transmission and reflection coefficients for high iv and fixed I follow from the 
asymptotics of the Bessel functions as in subsection 4.2: 



\R\ = I cos(n 7r)| 

|r| = |sin(n7r)| (5.12) 

The u dependence of the total cross section for transmission can be derived along the lines 
of subsection The leading terms in the potential (^^) at r ~ are 

(r+r_)'^-3 d-2 3 (r+r_)2(^ 
For fixed I the potential has a maximum at 

Tpeak - ■ (5.14) 

Equating y(rpeafc) with cj^ we get the transition value for I: 

Ic ~ . (5.15) 

Finally, we evaluate the total transmission cross section by employing qualitative classical- 
mechanics considerations. Classically the cross section is obtained by estimating the max- 
imal value of impact parameter b for which a significant transmission (or scattering, in the 
more general context) still occurs. Using b = l/uj we get 6max ~ Ic/^- Since in d dimensions 
the cross section scales as b'^~^, we find 

^~&max~ (^^j ■ (5.16) 

For Schwarzschild black holes the leading behavior near the singularity (at r = 0) is 



id - 2) 
1 



d-3 



/eff 



Thus independently of dimension one has c = 1/4 and hence total reflection: 



(5.17) 



n = 

T = . (5.18) 

5.2 The 2d black hole 

Black holes in two dimensions are apparently outside the pattern of the d > 4 black holes 
discussed in the previous sections (although as we shall recall Q they turn out to share 
similar features). Such solutions require, for instance, the presence of a dilaton. In this 
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subsection we shall consider two dimensional dilaton gravity with an Abelian gauge field 
which is inspired by string theory. The action is ^'^ 



(5.19) 



where g^y is the (string frame) metric, = 0, 1, $ is the dilaton, F^^ is the field strength 
of an Abelian gauge field A^, and A is the cosmological constant. In this two dimensional 
theory, the charged black hole solution is given by |3C] 

A(x -—dt 
r 

$(/>) = $0 - , (5.20) 

where ^>o is a constant and 

f{r) = 1 + ^ , r = e" . (5.21) 

Note that here ds^ oc ^ — / dt^, which is different from the generic d > 4 cases in eq. 
( |5.1| ). Nevertheless, the geometry of this 2d black-hole is similar to a two dimensional slice 
of the Reissner-Nordstrom solution, whose complete Penrose diagram is given in figure |2|. 
In the 4d case every point in figure ^ is actually a two sphere, while in the 2d case there is 
a non-trivial dilaton instead. 

An uncharged scalar field is minimally coupled to the background above as follows: 

S = J d^x^/^e-^* (^R + Ag'^'df.^dy^ _ If^ - A - gf"" df^ipd^iP - m^^) , (5.22) 

where in this section tp denotes the scalar field to distinguish it from the dilaton <I>. The 
wave equation of ip is: 

2$ 

V = -T= e-^-f g^^^^ a, V = {drr^fdr - r'dtt) V = m V , (5.23) 

where we used grr = l/(/?'^) and \/—g = e^^*~*°^ = 1/r. Hence, in the two dimensional 

case the dilaton $ plays again the role of the 4d spherical coordinates, now in the wave 
equation. 

As in the d > 4 cases, we separate variables 

V'(r,i) =^^0.(06*"* . (5.24) 

The equation for ipi^ is: 

{Oy/dr + uj^f-^ -m^)^^ = . (5.25) 



^^This can be obtained, for instance, as (part of) a low-energy effective action of an heterotic string in 
two dimensions (for a review, see Q). 
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Again, we define a coordinate r*, 



and a new field variable 



dr* 



dr 



rf{r) 
ip* := -y/r^ 



such that eq. ( 5.25 ) turns into: 
The effective potential is now 
where 

Next we inspect the behavior of the wave equation (|5.25| ) near the singularity of the 2d 
RN-hke black hole. When Q / 0: 



(5.26) 
(5.27) 

(5.28) 

(5.29) 
(5.30) 



(1 + + 0{r^)dr -% + {! + ...)^u;2 



(5.31) 



where stand for higher orders in r. As in the 4d case (^.6|), this equation is regular at 
r = for every uj. Moreover, near r ~ we have 



2Q2 



and hence 



(5.32) 



(5.33) 



16 r* 2 

On the other hand, for Q = - the Schwarzschild-like 2d black hole ~ the equation for 
is: 



-^{drr^fdr + u;^f-'-m^)i;^ 



(5.34) 



Again, as in the 4d case ( |3.9D , it is not regular at the singularity, but instead it is regular- 
singular. In the uncharged case, to leading order near r ~ we have 



and hence 



2M 



1 



(5.35) 



(5.36) 



As discussed in previous sections, the solutions to high frequency scattering wave equations 
(5.28) are given in terms of Bessel functions, which depend on the effective potential. In 
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the charged case we read the order of the relevant Bessel function from ( 5.33 ) and appendix 



n| = i , (5.37) 



and hence at high energies we find ( ^.12| ): 



irl = \R\ = ^ . (5.38) 



On the other hand, for the uncharged case, from (Ipl) and (|1|) we read: 



n = 

\R\ = 1 , T = . (5.39) 

To summarize, we see that the features of the 2d black holes discussed in this subsection 
are very similar to the 4d case. 

The 2d background ( ^.20| ) is also an exact Conformal Field Theory (CFT) background 
in string theory; in it was obtained from a family of ^^C^'^^^^^i^) quotient CFT sigma 
models by a Kaluza-Klein reduction^^. In string theory one is forced to include the regions 
beyond the singularity |3^ (this is argued, for instance, by using T-duality; for a review, see 
p^). The structure of the parent 5L(2,M) group allows one Q to find the exact solutions 
to the wave equation^^. In particular, it gives the exact reflection coefficient 



cosh I — - 1 cosh(u;) 

= / X ' (5-40) 

cosh ' — - — 



where a; is proportional to the energy of a massless particle scattered from the region beyond 
the black hole singularity, and r± are the locations of the event and Cauchy horizons: 

r± = M± - Q2 _ (5 41) 

Indeed, for the charged case — > ^ in the limit a; — > oo, and for the uncharged 
case l-Rp = 1 for every energy. These results are in agreement with ( 5.3S ) and ( 5.3S| ), 
respectively. 

6. Adding perturbations 

Until now we assumed that (j) was a scalar field with no mass, no charge, no back-reaction 
and no other interactions beyond minimal coupling to gravity. Now we will check whether 
our result that the time evolution is well-defined is disrupted by any of these perturbations. 
To our surprise we do not encounter any problem. 



^^In the bosonic case it is a semi-classical approximation, while in the superconformal extension this 
background is claimed to be exact 

^^Scattering waves are determined from vertex operators in the SL{2, R) CFT. Those are given by matrix 
elements in a unitary representation of the group. Equivalently, these are solutions to the Laplace equation 
(5.23), which in this case is a solvable hypergeometric equation. 
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6.1 Mass and charge 

Adding a mass m and charge q for the field (j) is completely captured by the form of the 



geodesic potential (2.19), namely one has 



T/wave f^geodesic , at^ 



where y^^'^^^^'^ [g simply y?^°'^'^^^'^ from ( 2.19| ) with the substitutions P ^ 1(1 + d — 3) and 



E ^ Lo, and AT4ff is the same as in the zero-mass zero-charge case ( 3.16| , ^^ ). 



Our analysis relied on the leading terms near r ~ 0, which we will now find to dominate 
over the added perturbations. For Q these are i/r^'^-^^ from AVcs and /^/r^"'"^ from 
^^odesic^ Since both m and q contribute at order l/r^'^~^, they are irrelevant near the 
singularity (for d > 4). Moreover, one can confirm that the form of the subleading terms 
still conform with ( |5.4[) and hence there are no log pieces in the solutions. 

Note however that since for q one cannot solve for u;^ in the wave equation and 
view it as eigen- values of some operator, any operator approach and considerations of 
self-adjointness (see section |^ would need to be changed. 

Finally, for Q = we might as well take q = and the leading terms are from 
AT4ff and P/r'^^^ from y^"'^'^^^'^ ^ while m contributes at and is again irrelevant. 

6.2 Interactions 

We may add non-quadratic terms to the action, namely adding terms non-linear in (f> to the 
wave equation. Much of our previous analysis, especially the separation of variables relied 
on linearity, hence we should reconsider it, starting from the non-linear wave equation 

= a(P-gU{ct)) , (6.2) 

where g is a coupling constant and U{(p) is a non-linear potential term. One can attempt 
to solve this equation by perturbation theory with g being the small parameter, namely 
one expands 

(l) = J2g^ (j)^^^ . (6.3) 
j 

We start at zeroth order with a solution cj)'^^^ to the linearized equation 



□(/)(°) = (6.4) 
and consider whether the first correction is regular at r = 0. The first correction satisfies 

=gC/(</,(0)) . (6.5) 
In order to solve this equation we separate the angular and time variables 

r^n^ict^^^} =gr^ \U{<P^H , ■ (6.6) 
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bmce 

^(0) 

is regular by assumption, so are the components This is a non- 

homogeneous ordinary (second order) differential equation. A solution in the vicinity of 
the singularity may be obtained by Taylor expanding the equation (see also section |^. 
The homogeneous solutions are the smooth (fy^^j , while the non-homogenous term added to 
( |5.3| ) is O {f~^{r)^ and hence it is subleading and does not change our ability to obtain a 
(univalued) solution. This will happen at higher orders in g as well. 

6.3 Back-reaction for RN in d > 4 

So far we worked in the linear approximation where one neglects the back-reaction of 
the scalar field on the background geometry. Let us check whether including the back- 
reaction disrupts our result, namely whether we can find any obstruction for the solutions 
to the linear equation from being extended to solutions for the full non-linear system of 
background plus scalar field. Our analysis will be limited to lowest order back-reaction, 
and we shall not cover all cases, but the indications are that these solutions do survive 
back-reaction. 

Let us define a small parameter e such that cp is first order in e. Then at second order 
the background may change - the metric due to a T^u source, and the electromagnetic (EM) 
field due to the charge associated with cj). For simplicity we may consider cp to be neutral. 
Since the metric is singular to start with, we cannot use the straight-forward criterion 
that the background remains regular after back-reaction. A conservative view would be to 
continue and look at third order whether the changes in as a result of the changes in the 
background make it singular. Another approach is to require that the singularity in the 
metric does not become worse after back-reaction. The latter has the advantage that one 
can stop at second order. 

Moreover, for a charged black hole we can even save us the work of computing the 
second order by the following observation. Let us compare the stress-energy tensor of the 
scalar field with that of the EM field, namely the background source. The action for gravity 
-|- EM -|- an uncharged scalar field is 

S = j (fx^g {r - -^F^ - g^''d^<pd,(l) - m^cj)''^ , (6.7) 

where g^^ is the metric, F the field strength and 4> an uncharged scalar field of mass m. 
The stress tensors of the gauge and scalar fields are 

= d^4>du4> - \g^u [{d^? + m'ct?] (6.8) 
We start by comparing the scalar quantities 

{d<t>? = dr<t>dr^g--^^ 

{F') = KtFr,g-g''^-^^ (6.9) 
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where we use the background ( p.lD and the property that dr4> are regular at the singu- 
larity. We see that close to the singularity the scalar field seems to add a negligible source 
on top of the one from the EM field. 

We now proceed to compare all components 

We see that indeed for the rr, tt and 6*0 components T^i^'* ~ T^f ^ is negligible, while T^f^ 
is regular. Altogether we consider this to be strong evidence that back-reaction is weak 
and does not change the smoothness of the solutions. 

For a charged field (j) the simple argument above does not seem to work. In principle 
one should determine the back-reaction to the metric and gauge field and whether they are 
more or less singular than the original background. A preliminary analysis turned out to 
be involved, and so we did not reach any conclusions for this case. For similar reasons we 
do not discuss the back-reaction to the negative-mass Schwarzschild either. 

6.4 Back-reaction for the 2d black-hole 

In two dimensions, the action for an uncharged scalar field coupled to the metric and 



dilaton is (see subsection 5.2): 



S = j (fx^e'^"^ (^R + 4g'"'d^<^d^^ - ^F^ - X - gf""" d^tpd^i; - m^^-^^ , (6.11) 

where gfj_iy is the string frame metric, $ the dilaton, F is the field strength of an Abelian 
gauge field A, i/j is the scalar field (to distinguish it from the dilaton) and A is the cosmo- 
logical constant. 

Using the background ( ^.20 ) and the property that ip, drijj are regular at the singularity, 



we find the following behavior of the scalar stress tensor ( |6.8D near the singularity 

e^'^T^P ~ 1 



The dilaton stress tensor is 



e2*r^*) = 4V^V,$ - 45^, (v2$ - {d^f - Ia) (6.13) 



and has the following leading behavior near the singularity 



,2<I.y(<I.) _ ]_ 



2*ri*) = (6.14) 
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Again, there is no strong back-reaction to the metric, since component- wise each component 



of Tj^tt"^ is either subleading to Tj^u' or regular. 

To check the back-reaction to the dilaton, we define 



X 



(6.15) 



in terms of which the action becomes 



/ 



Sx\/ —g 



(6.16) 



The equation of motion for x is 



1 



AUx - [R- -F' -\- g^'d^ijd.i, - ) X = , 



(6.17) 



where □ is the Laplacian without the dilaton. In principle, we should expand the so- 
lution to ( |6.17] ) to second order in the perturbation and compare the correction to its 
initial value However, just like for the metric, we recognize that the new source for x, 
[g'^'^ d^ipdu'ip — m?"ilP), is negligible at the singularity relative to the background sources 
{R — \F'^). Therefore, we do not expect a strong back-reaction for the dilaton either. 
Summarizing, we find indications for a weak back-reaction near the singularity of the 2d 
charged black-hole. 



7. General conditions for resolution of a singularity 

7.1 Uniqueness of decomposition and natural boundary conditions 

In order for the "singularity smoothing" procedure which we defined in section ^ to make 
sense we need the existence and uniqueness of decomposition of any initial condition into 
a linear combination of eigen- functions (j)^^. We know that Hermitian operators have this 
property, namely their eigen-functions form a basis of function space. Thus we are inter- 
ested in representing the radial equation U^^i cpui = ^ through an operator and studying its 
Hermiticity, and so we define the operator L by 

f\J^i = u? - L[dr,r] 

L:=-f,drfr'dr + ^^±^ (7.1) 

and we specialize back to 4d for definiteness. We note that L is not regular at r = 
even for RN but only due to an overall factor of relative to the regular radial equation 
f~^\I\^i (p.5|). One may wonder whether the decomposition property would nevertheless 
hold for L. This will be answered in the negative as we shall show that the functions (f)^ 
do not span all possible functions. 



17 / 



At zero order we have x = XoxA" (5-20, 5.21). The second solution to (6.17) is subleading near r = 



so turning it on at higher orders would not change the result. 
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Since the differential equation (|3.5| ) is smooth at r = one can Taylor expand the 
equation and solutions 

oo 

cPUr) = ^4>^fr^ . (7.2) 

3=0 

Substituting this in the differential equation transforms it into a recurrence equation for 
the Taylor coefficients (pul . 

For RN one notes that the coefficient of is which is O(r^) and hence the 
recurrence equation is w-independent up to order O(r^) (inclusive) where (p'^^ is determined. 
Among the first 6 Taylor coefficients cj)!^^ , • • • , 4>u^ the first two may be considered to be 
initial conditions for the recurrence equation and the other four may be considered to 
be a function of them. In other words we find 4 ^-independent linear relations among 
\ ■ ■ ■ Owing to the w-independence of these relations they continue to hold for 

any function which can be expressed as a linear combination 

(pi = J dio ((^+(cj) +(/>_(cj)) (P^{r) , (7.3) 

where (f)±{uj) are arbitrary coefficients. This means that there are 4 conditions on the first 

6 coefficients of any function in span(0i^). 

For Schwarzschild the phenomenon is similar but the details are different. We have 
= O(r^) and hence the recurrence relation is free up to (and including) the 4th coeffi- 

(3) 

cient 0ij . Since we consider only the regular solutions then the recurrence relation requires 
a single initial condition cp^^ while (/>i^^ , • • • , i^af ^ may be considered to be an w-independent 
function of it, and hence there are necessarily 3 conditions on the first 4 Taylor coefficients 
of any function in span((/)(^). 

We see that span((/>i^) does not contain all functions, but rather there are constraints on 
the initial conditions at r = 0. This is why we required the time-evolution to be well-defined 
only for wave packets prepared far away, but otherwise arbitrary. We expect that once the 
domain of L is correctly defined it will be Hermitian and uniqueness of decomposition will 
hold and consequently the uniqueness of time evolution. The mathematical term that we 
expect to hold is that L is "essentially self-adjoint" namely that it has a unique Hermitian 



(i.e. self-adjoint) extension (see for a text book and fl^, 



Initial conditions at r=0 

Let us find some necessary boundary conditions at r = in response to the "non-span" 
property just discussed. Since L defined in ([7.lD is singular at r = du may be ill-defined 
(which reflects the singular nature of the wave equation ( |3.2D ). Since cj) is finite at r = 
we need dtt(t> to be finite as well and thus we clearly need 

L(/> = finite (7.4) 

where (p can be either of the initial conditions 0j, (pi. 



In the rest of this subsection "finite" should be understood to mean "finite at r = 0.' 
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For RN from the leading behavior for /rn (|2.2| ), we get that for a regular and generic 
(j) the Laurent expansion of Lcf) starts at degree (—4), and thus eq. ( |7.4| ) represents 4 
conditions (for degrees —4 < deg < —1). Moreover, if we Taylor expand 4> = Y2JLo 4'^''^ 
then at order deg the largest j for which (j)^^^ appears in the equation is j = deg + 6, and so 
we have 4 (linear, homogeneous) conditions for the 6 coefficients (f)^^^ , ■ ■ ■ , (t>^^^ ■ These are 
actually the conditions found in the paragraph around eq. ( |7.2| ), which limit our choice of 
initial conditions in the vicinity of r = 0. 

Similarly, for Schwarzschild (^^) represents 3 conditions (for degrees —3 < deg < — 1) 
among the first 4 Taylor coefficients which agree with the findings in the paragraph after 
eq. (P). 

However, there are additional requirements: the constraint ( [7.4[) must be compatible 
with the time evolution according to the wave equation 

= -du(l> (7.5) 

namely we must require 

att[L(/>] = finite (7.6) 
By the wave equation (^^) this is equivalent to 

= finite (7.7) 

This gives 4 additional conditions (in addition those of ([7.4|)) involving the first 12 Tay- 
lor coefficients for RN and 3 additional conditions for the first 8 Taylor coefficients for 
Schwarzschild. 



Now we may readily generalize, and find 
Necessary boundary condition: 



finite Vn > 1 



(7.8) 



(at r = 0). 



This condition holds both for RN and Schwarzschild and is explicitly compatible with 
time evolution. Note that it can also be interpreted as 9^ (/> = finite Vn > 1. It is plausible 
to us that this condition is not only necessary but also sufficient for a unique time evolution 
as well since showing that (p has finite time derivatives of any order at the initial moment 
of time comes close to providing a well-defined time evolution, but we shall not attempt to 
pursue this point. 



As we already mentioned, the boundary conditions (7.8) are automatically satisfied 



for characteristic b.c. (at past null infinity) and more generally b.c. which vanish at a 
neighborhood of r = 0. 

Another perspective is to consider a numerical computer implementation of this time 
evolution. One may either have no grid points at r = in which case the time evolution 
is well-defined, or put a grid point at r = and find its time increment as a limit of the 
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time increments of neighboring points. The singularity will tend to "expand" numerical 
errors, but those should be possible to tame by decreasing the grid spacing. It would be 
interesting to analyze this further and/or to perform the implementation and determine 
the behavior at r = 0. 

7.2 General conditions for wave-regularity 

Having shown that Reissner-Nordstrom is wave-regular allowing for transmission across 
the singularity, we would like to abstract the general conditions on a spacetime to have 
such a resolution. 

The crucial property of the RN singularity is that the two solutions (f)^^ have a unique 
continuation across r = 0, while a priori there could have been multi-valued functions such 
as the log's which appear in the Schwarzschild case. Hence a wave packet constructed 
from the (p^i has a prospect of crossing the singularity in a well-defined manner. This 
is true not only in 4d RN where the time-separated ODE is regular but also for d > 4 
RN where the ODE becomes regular-singular rather than smooth, but nevertheless the 
two characteristic exponents arc integral and there arc no log pieces so that the functions 
(j)^^ are univalued in the vicinity of the singularity. The condition on the characteristic 
exponents can be generalized to allow for a reparanictcrization of coordinate r ^ f ~ r'^i 
and a linear redefinition of the field (p ~^ "f"^^ 0^ for some constants ci , C2 . In general 
can be series expanded as = Z^fcLo 4>^'^^ where p = pi,2 is one of the two 
characteristic exponents and A is the "series step-size" (A = 1 when the functions in 
the differential equation are meromorphic) and cfP^'^ are some constants. The quantity 
{pi — /52)/A is invariant under the double transformation above, and hence there exists a 
transformation such that all three pi, p2 and A are integral and the function is univalued 
exactly if the ratio above is rational. 

We summarize the above by 
Necessary general condition: the "eigen-functions" (f)^ should be univalued in the 
vicinity of the singularity. This is equivalent to: 

• The equation for (f)^ is either regular or regular-singular. If it is regular-singular we 
also require the next items: 

• The difference of characteristic exponents is commensurate with the "series step-size" 
(defined in the previous paragraph) 

P^^^ G Q . (7.9) 

• There are no log pieces in the solutions. 
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8. Summary and discussion 

In this paper we saw that physical predictabihty can be restored in the neighborhood of 
the charged (RN) black-hole singularity as well as negative mass Schwarzschild, once one 
considers waves rather than particles. For RN this is done by gluing two spacetimes over 
the time- like singularity, thereby adding another asymptotic region. An observer at infinity 
in the additional region views a spacetime with a negative mass. Several alternatives exist 
for the total Penrose diagram depending on the relative size of and Q^. This singularity 
appears to have the physical nature of a beam splitter. 

For negative-mass Schwarzschild there is a natural "regularity" boundary condition at 
the singularity cutting off spacetime there, and thus it can be physically interpreted as a 
perfectly reflecting mirror. 

We subjected this picture to several tests by adding perturbations in section ^. Con- 
sidering a field with both mass and charge we found that such terms are subleading at 
the neighborhood of the singularity. Considering added interactions the field equation be- 
comes non-linear and we were satisfied in confirming that there is no obstruction (such as 
creating a divergence) in extending the linear solutions to next order. Similarly we then 
considered the effect of non-linearities from back-reaction, and found that in some cases a 
simple argument protects the regularity of (f) at the next order. Altogether the results were 
surprisingly resilient to perturbations. 

It would be interesting to perform some additional tests and generalizations: 

• Study fields with other spin. Of particular interest are the electromagnetic and 
gravitational fields, and also the Dirac ,3 = 1/2 field. A preliminary investigation of 
the electromagnetic and gravitational fields in d=4 indicates that their polar modes 
behave just as scalar fields, i.e. 75% refiection and 25% transmission for high energy 
at fixed /. The axial modes appear to be more subtle, however. 

• Study other backgrounds such as extreme RN and rotating black-holes. 

Rotating black-holes deserve a special discussion here. They may be linked to our 
work in two different ways. First, the motivation to the present analysis partly emerges 
from the desire to explore the physical phenomena that may take place deep inside realistic 
rotating black holes. From this point of view the spherical charged black hole serves as 
a toy model for the more complicated, non-spherical, spinning black hole. In the second 
link we can employ the spinning black holes to test the singularity-resolution approach 
developed here. In the Kerr-Newman solution the r = singularity forms a ring rather 
than a hypersurface. Once the field is decomposed into spheroidal-harmonic modes, the 
field equation for each mode is perfectly regular at r = 0, hence there is no doubt about the 
proper continuation. Now, when the spin parameter a is taken to vanish, the ring's radius 
shrinks to zero, and the spacetime becomes RN. One may therefore define the RN extension 
of the field beyond the r = singularity to be the limit a ^ of the corresponding field 
in the Kerr-Newman case. The obvious question is, therefore: does this procedure yield 
exactly the same extension as that constructed above? At least for a scalar field in 4d RN 
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the answer is found to be positive. It still remains to check whether this is also the case 
in the Q = case. Namely, in the uncharged Kerr case, when a — > and the spacetime 
becomes Schwarzschild, does one recover the full-reflection b.c. advocated above? This 
still needs to be verified. 

We would like to mention several other issues at the classical level: 

• In all of our examples one of the spacetimes has a negative mass and a globally- 
naked singularity. Such physical objects raise problematic issues, and we name only 
a few: anti-gravity, acceleration reversed to force, and also inconsistency of their 
construction with the Cosmic Censorship conjecture (see "Perspectives on r < 0" in 
section ^). 

It would be interesting to determine whether Nature allows the actual construction 
of any spacetime with a wave-regular timelike singularity. The mechanism that im- 
mediately suggests itself would be to reconsider a gravitational collapse of a charged 
spherical shell (in the positive-mass universe), leading to a RN black-hole geometry 
(outside the collapsing object). However, the instability of the inner horizon raises 
some doubts about whether the RN-like r = singularity will indeed form in this 
process. 

• Implications for Cosmic Censorship. At the very least it shakes its rationale since 
given "reasonable" initial conditions Cosmic Censorship is supposed to "protect us" 
from loss of predictability, namely of losing unique time evolution, by forbidding 
naked singularities, but here we see that such singularities may not mean the loss of 
predictability after all. 

• Since for RN we advocate a picture where there are two spacetimes which are consis- 
tently glued at the singularity — and since a macroscopic measuring device cannot 
cross the singularity (at best it will bounce back, but it may also be destroyed by 
tidal forces, or be "beam-splitted," in the worse case) one may take two different 
points of view on the physics corresponding to either of the two observers on the two 
sides of the singularity. Namely, the story of the resolution of each singularity will 
be told in two different versions corresponding to the two observers (this applies to 
both types of RN spacetime, namely \Q\ < \M\ and \Q\ > \M\ - see figures §,||). 

Going beyond classical GR there are open questions as to the quantum gravity consis- 
tency and properties: 

• We must note that our wave equation is outside its domain of validity near the 
singularity due to the presence of high curvature. However, a preliminary analysis 
shows that quantum corrections do not produce essential singularities in the equations 
and start altering the field only at a Planck (proper) distance from the singularity. 

• Semi-classical quantization and Hawking evaporation (under study). 
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• As one passes to the quantum theory one may suspect that the pathologies associated 
with negative-mass spaces to only grow worse. Therefore, we need to be very cautious 
in discussing possible implications for "legitimizing" spacetimes with such singularities 
as possible solitons, namely which spacetimes should be considered to contribute to 
the path integral as admissible saddle points. Here we took an open-minded approach 
of exploration and we hope that the various interesting issues which get raised will 
be studied further. 
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A. Kruskal-Szekeres coordinates 

For completeness, we recall here the derivation of the metrics in the Kruskal-Szekeres 
coordinates and the definition of the functions g{r) which are mentioned in the text and 
are used to implicitly define r in these coordinates. 
First one defines the "tortoise" coordinate r* by 

dr* := J . (A.l) 

r* diverges on horizons but is finite at r = 0. Then one passes to light-cone coordinates 

V ■= t + r* 

u:=t-r* (A.2) 

and finally in the neighborhood of a horizon with surface gravity k (for Reissner-Nordstrom 
we have k± = (r+ — r_)/(2r^) while for Schwarzschild k = l/(2ro)) the Kruskal-Szekeres 
coordinates are given by 

V := ibexp(ibKf) 

U := ±exp{±Ku) (A.3) 

and the sign inside the exponent is chosen such that the horizon is located at UV = 0, 
namely it is X-shapcd. 

In these coordinates the metric is given by 

ds^ = — ^ dU dV + dft^ , (A.4) 
K^g 
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where the function g{r) is defined by 



g{r) := exjp{±2Kr*) (A.5) 



and the sign is chosen to be the same as that for V in ( [A.3| ). r is not a coordinate anymore, 
but rather it is imphcitly defined in terms of the U, V coordinates by 

-UV = g{r). (A.6) 

Note that by construction g{r) has a zero at the horizon and thus the horizon is exphcitly 
smooth in these coordinates since in the prefactor of dU dV the zero in / gets cancehed 
against the zero in g. 

In Reissner- Nordstrom there are two horizons at r±, and accordingly two Kruskal- 
Szekeres-hke planes for {U±, and two functions g±{r). The "upper plus" quadrant 
(C/+, y+ > 0) is glued to the "lower minus" quadrant ([/_, V- < 0) according to 

[7+/^+ = i-U-)^/"- (A.7) 

and similarly one may continue and glue the upper minus quadrant to another copy of the 
plus plane, creating the maximal analytic extension of Reissner-Nordstrom which consists 
of an infinite chain of alternating plus and minus planes. g±{r) are given by (see figure |l]) 

5_(r) = (l-—) (l-—) exp(-2/t_r) (A.8) 



r_ / \ r_|_ 
While for Schwarzschild we have (see figure |^ 

gschw{r) = ( — - 1 ) exp(r/ro) . (A.9) 

Finally, in order to get the Penrose diagram (see figures one customarily takes the 

conformal transformation 

Up = tg-\U) 

Vp = tg~\V) (A.IO) 
although any other transformation which maps the real line to an interval is admissible. 

B. Regular singularities of ordinary differential equations 

Let us briefly recall the definitions. A linear second order differential equation 

[a{r)drr + b{r)dr + c{r)](j) = (B.l) 
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is regular-singular at a point ro, which we will assume without loss of generality to be 
ro = 0, if after normalization such that a(0) = 1 there are singularities at r = in b{r) or 
c(r) of limited type: b{r) may have at most a first order pole, and c(r) at most a second 
order pole. This guarantees that the solutions will have at most poles or branch cuts at 
r = but not an essential singularity. More specifically the two solutions have the leading 
behavior 

<^i,2 ~ rP''' , (B.2) 

where pi^2 are called the characteristic exponents, and they are the solutions to the 
quadratic equation 

a{ro)p (p - 1) + 6(ro) rp + c(ro) = (B.3) 

gotten from substituting the leading behavior above into the equation. When pi = P2 = P 
the leading behavior is (pi = r^, (j)2 = log(r). Finally, one can expand the solutions into 
a series tj) = rP X^j^o ^^^^ where (j)'^^'^ are some constants, except that in the case when 
P2 — pi is a positive integer </>i may contain also a piece proportional to log(r) (j)2 and thus 
contains a log. 



C. Basis functions for the characteristic formulation 

Here we define the basis functions , ^ -j for the characteristic formulation (see 

subsection ^7^ ). Since the radial equation is a second-order ODE, the radial functions 
(r*) for a given a; form a two-parameter family. Owing to the asymptotic behavior 
of the radial functions, Eq. (3.23), we may choose two basis functions <;^^io)('^*) ^"^^ 
(/)*^Q j^j(r*), defined as follows: (/"^(i o)('"*) is the radial function which at r* — > +oo (the 
inner horizon) has the asymptotic form 

C{i,o) - r(a;)exp(-iu;r*) 

and at r* ^ — oo (the negative-r asymptotically- flat region) has the asymptotic form 

C(i,o) - exp(-i ur*)+ i?(a;)exp(+i u; r*), 

where T{uj) and R{ijo) are two unconstrained coefficients (these coefficients turn out to be 
the transmission and reflection coefficient; see section Q). Similarly, 4>Ij^q i)(^*) the radial 
function which at r* ^ — oo has the asymptotic form 

C(o,i) - r'(a;)exp(+iu;r*) 

and at r* — > +oo has the asymptotic form 

C(o,i) - exp(+i u r*) + R'{uj)exp{-i uj r*), 

where again T'{uj) and R'{uj) are two unconstrained coefficients. Now, the two desired 
solutions (j)*^^ and (p*^^ ^ ^ are simply given by 

'/'K,o)(^*'*) = C„(i,o)(^*)exp(«w„t), 
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As one can easily verify, these two solutions take the asymptotic forms 
^K,o) - T{uj^)exp{iuj^ v) (r* ^ +00), 

'^lujv 0) ~ exp(ia;^ v) + R{io)exp(i lOv u) (r* — > —00), 

and 

^(0,c^„) - T'{uJu)eyiv{ii^uU) (r* -00), 

^lo.i^u) - ^^P(^ + -R'(^^)exp(i v) {r* +00), 

in agreement with the above definitions of 61, , n^ and (hin , , v 

D. Bessel functions 

Let us assemble a few useful properties of the Bessel functions. A Bessel function of order 
n satisfies the equation [dxx + i^/^) dx + ^ — n'^/x'^] Jn = 0. It will be more convenient for 
us to use the equivalent definition 



dxx + I 1 



- 1/4 



x"^ 



V^Jnix) = . (D.l) 



For small x 

Ux) = x^ Aiix^) , (D.2) 

where A;((x^) are certain functions of analytic at x^ = 0. The asymptotic behavior is 
better described by 

= (Jn cos(n7r) — J-n)/ sin(n7r) (D.3) 

and one has 

J„ ± iy„ ~ v^2/(7rx) exp [±i{x - n7r/2 - 7r/4)] (D.4) 

as X — >■ +00. 
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